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Abstract 

We study the role of density-dependent low-momentum nucleon-nucleon interactions 
in describing the anomalously long beta decay lifetime of C. We approach this problem 
both from the perspective of chiral effective field theory, in which genuine three-body 
forces generate an effective density-dependent two-body interaction, as well as from the 
perspective of Brown-Rho scaling, in which the masses and form factor cutoffs in one- 
boson-exchange interactions are modified in a dense nuclear medium due to the partial 
restoration of chiral symmetry. The beta decay transition of 14 C to the ground state of 
N is calculated within the shell model using a model space consisting of two Op-shell 
holes within a closed 16 core. The effective Op-shell interaction is calculated up to second 
order in perturbation theory with single-particle energies extracted from experiment. We 
find that both three-nucleon forces and Brown-Rho scaling medium modifications give 
qualitatively similar results not only for the ground state to ground state Gamow- Teller 
transition but also for Gamow- Teller transitions from excited states of 14 C to the ground 
state of 14 N. In this way, it is observed that at a low-momentum scale of Ai ow _k = 2.1 
fm _1 , medium- modifications of the nuclear force play an essential role in increasing the 
lifetime of 14 C from a few minutes to an archaeologically long one of 5730 years. 

PACS: 21.30.Fe, 21.60.Cs, 23.40.-s 
Keywords: Effective field theory at finite density, chiral three-nucleon force, Brown-Rho scaling, shell- 
model calculation of 14 C beta-decay. 

1 Introduction 

In recent years our understanding of nuclear structure physics, and the nuclear force in par- 
ticular, has benefited greatly by incorporating aspects of the fundamental theory of strong 
interactions, QCD. Although QCD is highly nonperturbative at the energy scales characteristic 
for nucleons bound within nuclei, one can exploit the symmetry breaking pattern of QCD, 
namely the explicit and spontaneous breaking of chiral symmetry, to constrain the underly- 
ing effective theory for interacting nucleons. At low energies, this is accomplished through 
chiral effective field theory [1]-[10], in which the light pseudoscalar Nambu-Goldstone bosons 
(e.g. pions) generated by the spontaneous breaking of chiral symmetry comprise the low-energy 
mesonic degrees of freedom relevant for nuclear interactions. All other heavier mesons are 
integrated out and their effects are absorbed in a series of short-range contact interactions. 

The above description utilizes the spontaneous breaking of chiral symmetry at low energies. 
However, as the temperature and/or density of a hadronic medium increases, the order param- 
eter for chiral symmetry breaking (the scalar quark condensate (qq)) decreases in magnitude 
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and is expected to vanish at an energy density around 1 GeV/fm 3 , signaling the restoration 
of chiral symmetry. The vacuum structure of QCD, encoded in various quark and gluon con- 
densates, is expected to play an important role in determining the properties (e.g. masses) of 
light hadrons [1 1]— [13] , with the scalar quark condensate being of primary importance. The 
work of Brown and Rho [12, 14] was one of the first attempts to connect medium-dependent 
meson masses to the partial restoration of chiral symmetry as the temperature and/or density 
changes. By imposing scale invariance - a symmetry of classical, massless QCD - on chiral 
effective Lagrangians, Brown and Rho obtained a scaling law in which the masses of hadrons, 
except the pion, decrease with the cube root of the scalar quark condensate. 

A natural testing ground for Brown-Rho scaling is heavy ion collisions, which produce 
hot/dense hadronic matter in the vicinity of chiral restoration. For a recent review see refs. 
[15, 16]. However, even for densities at the center of heavy nuclei, the scalar quark condensate 
is expected to decrease in magnitude by about 30% as found from model studies [17], low- 
density expansions [18], and calculations within in-medium chiral effective field theory [19]. The 
decreasing of meson masses at these densities would have important consequences for nuclear 
structure, since the intermediate-range attraction and strong short-distance repulsion in one- 
boson-exchange models of the nucleon-nucleon potential arise largely from the exchange of the 
a and oj mesons, respectively. Moreover, phenomena that are sensitive to the tensor component 
of the nucleon-nucleon interaction would be affected by the decreasing of the p meson mass. 
It was therefore hoped that evidence for the partial restoration of chiral symmetry at finite 
density could already be found by studying the properties of normal nuclei. 

Brown and Rho proposed employing the hidden local symmetry approach [20, 21] as a 
framework for studying mass scaling effects. In this formalism vector mesons, such as the p and 
oj, are treated as gauge bosons that acquire mass through the Higgs mechanism. Presently a 
consistent hidden local symmetry approach incorporating baryons has yet to be fully developed. 
This limits the applicability of the formalism to dense nuclear systems. A more practical 
approach has been to directly modify the masses of mesons in one-boson exchange interactions 
[22]-[26]. This leads to a density-dependent nucleon-nucleon interaction that can be used in 
nuclear structure studies. One aims to achieve a qualitative understanding of the effects of 
in-medium hadron masses on nuclear structure. 

Not only can the mass of a hadron be affected by a background medium, but more generally 
the entire spectral function will be modified due to Pauli blocking effects, collisional broadening, 
and the opening of additional decay channels in a hadronic medium [27]- [31]. Indeed the study 
of spectral functions has become the dominant focus both experimentally and theoretically 
for probing the tendency towards chiral symmetry restoration below the transition point. The 
importance of resonance-hole excitations in modifying the spectral features of the vector mesons 
p and oj has been studied in detail. 

Chiral effective field theory, representing the interface between QCD and nuclear physics at 
low energy and momentum scales, takes a different approach by integrating out mesons much 
heavier than pions. The short-distance dynamics is encoded in contact terms. Intermediate- 
range effects that would emerge from a exchange and the p-exchange tensor force, are now 
generated by explicit two-pion exchange processes, subject to Pauli principle corrections in the 
nuclear medium. Resonance-hole effects are incorporated through genuine three-nucleon forces. 

With the development of highly-accurate methods for solving the nuclear many-body prob- 
lem for light nuclei, as well as infinitely-extended nuclear matter, it is now well established that 
free-space two-nucleon interactions alone are insufficient to describe the properties of dense nu- 
clear systems. In particular, calculations using only nonrelativistic two-body forces have been 
unable to accurately reproduce the binding energies and spectra of light nuclei [32]- [34], the 
saturation binding energy and density of isospin-symmetric nuclear matter [35]-[39], and finally 
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the nucleon-deuteron scattering differential cross sections at intermediate energies [40]-[43]. In 
all of these cases, the errors arising from the solution of the nuclear many-body problem are 
small compared to the uncertainties in the underlying nuclear force model, which suggests the 
need to include explicitly degrees of freedom beyond nucleons and mesons, or equivalently the 
need to include three-body forces. 

In this work we study the unique problem of the beta decay lifetime of 14 C, which has been 
shown [25, 44] to be highly sensitive to density- dependent effects in the nuclear interaction. 
Describing the beta decay lifetime of 14 C has long been a challenge in nuclear structure physics, 
because the anomalously small Gamow- Teller (GT) transition matrix element M GT ^ 0.002 
results from the cancellation between terms that are otherwise on the order of 1. The valence 
particles of 14 C inhabit a region with a large nuclear density, and therefore medium modifications 
to the nucleon-nucleon interaction could be sizable. In ref. [25] it was suggested that a decreasing 
in-medium tensor force from Brown-Rho scaling plays a significant role in the suppression of 
the Gamow- Teller transition. This was inspired by early work [45, 46] that highlighted the 
sensitivity of the beta decay lifetime to the tensor component of the nuclear interaction. In 
contrast, a later study including three-nucleon forces [44] found qualitatively similar effects 
that resulted from additional short-distance repulsion coming from the contact three-body force. 
Within Brown-Rho-scaled nucleon-nucleon interaction models such an effect is generated by the 
decreasing of the u meson mass. We wish to revisit the problem of the 14 C beta decay lifetime 
and suggest that, indeed, in both models of the low- momentum density- dependent nucleon- 
nucleon interaction, additional short-distance repulsion generates significant suppression of the 
Gamow- Teller transition. 

This paper is organized as follows. In Section 2 we describe the framework of chiral effective 
field theory, and how a density-dependent two-nucleon force can be constructed from the chiral 
three-body interaction by summing over the filled Fermi sea of nucleons. In Section 3 we 
motivate Brown-Rho scaling and study the effect that such hadronic medium modifications 
have on the nuclear force. In Section 4 we apply both models to the problem of the extremely 
long beta decay lifetime of 14 C, which was suggested to be particularly sensitive to the density 
dependence of the nuclear force. We highlight the important role played by the short- distance 
repulsion (increasing with the nuclear density) in both approaches to the in-medium nuclear 
interaction. We end with a summary and conclusions. 

2 Density-dependent nucleon-nucleon interaction from 
chiral effective field theory 

Chiral effective field theory has proven to be a highly useful framework for studying the nuclear 
interaction problem (for current reviews, see refs. [10, 47]). By implementing the spontaneous 
and explicit breaking of chiral symmetry in QCD, one constructs a low-energy effective field 
theory of interacting pions and nucleons order- by-order in powers of Q/A x , where Q represents 
a nucleon momentum or the pion mass. The long- and intermediate-range parts of the nuclear 
force are generated by lir, In (and more) exchange processes, whereas the short-distance dy- 
namics due to heavy mesons and baryon resonances are integrated out and their effects are 
encoded in nucleon-nucleon contact terms. Currently the nuclear interaction has been calcu- 
lated up to 4th order (N 3 LO) in the chiral power counting [6]-[9]. It is at this order (and 
beyond) that it is possible to achieve an agreement with empirical NN scattering phase shifts 
that is comparable to previous high-precision NN potentials [48]-[50]. By adjusting the 29 pa- 
rameters (mostly low-energy constants) that occur at this order, the 1999 database for np and 
pp elastic scattering up to E\ ah = 290 MeV can be fit with a x 2 /dof of 1.1 for np scattering and 
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1.5 for pp scattering. Furthermore, the experimental deuteron binding energy, charge radius, 
and electric quadrupole moment are very well reproduced by the chiral N 3 LO potential [9]. 

One decisive practical advantage of the chiral effective field theory approach is that two- and 
many-body forces are generated consistently within the same framework. Chiral three-nucleon 
forces first appear at N 2 LO. There are three different components, which we show in diagrams 
(a), (b), and (c) in Fig. 1. The long-range, two-pion exchange interaction is given by 

3 " " >h 8 # (* a + <m 2 + o 3 3 ' 

where qa = 1.27, /„- = 92.4 MeV, m n = 138.04 MeV is the average pion mass, cfi = p( — pi is 
the difference between the final and initial momentum of nucleon i, and 

Fg = 5 a ? (-4dm* + 2c3« • q 3 ) + c 4 e a ^T^ k ■ (q, x q 3 ) . (2) 

The three low-energy constants of the long-range 3NF, namely C\ = — 0.76 GeV -1 , c 3 = 
— 4.78GeV _1 , and c 4 = 3.96 GeV -1 appear in the two-pion exchange component of the NN 
interaction and therefore can be determined from fits to low-energy NN phase shifts [51]. The 
medium-range, one-pion exchange three-nucleon interaction is proportional to the low-energy 
constant cd/A x and has the analytic form 

v (i*)__ V 9ACp & r qj ( ) 

with A x = 700 MeV. The filled black square in diagram (b) of Fig. 1 symbolizes the correspond- 
ing two-nucleon one-pion contact interaction. Finally, the short-distance contact three-body 
force is proportional to the low-energy constant c^/A x and is given by 



V^P = E ■ r s . (4) 



In eqs. (l)-(4) the nucleon labels k can take the values 1,2,3 which results in six possible 
permutations for each sum. The low-energy constants of the medium-range (V 3 ^) and short- 
range (vffi) three-body forces are fit to properties of few- nucleon systems, such as the triton 
binding energy together with the nd doublet scattering length [52], the 4 He binding energy 
[53], or the binding energies and spectra of light nuclei [34]. The N 3 LO three-body interaction 
contains numerous one-loop diagrams but no additional low-energy constants and is currently 
under construction [47]. 

In this work we utilize renormalization group techniques to decimate the chiral N 3 LO two- 
body interaction down to an energy scale in the vicinity of A = 2.0 fm" 1 , following the details 
provided in refs. [54, 55]. Such a decimation procedure has the desirable feature that the 
resulting "low-momentum interaction", referred to as Vi ow _k, is nearly unique regardless of 
which underlying high-precision NN potential is employed, and therefore Vi ow -k provides a 
model-independent description of the nucleon-nucleon interaction below A = 2.0 fm _1 . The 
analogous evolution of realistic three-body forces down to this scale has yet to be achieved. 
However, since such a decimation would change primarily the short-distance part of the three- 
body interaction, a practical alternative has been to treat the low-energy constants and ce 
as functions of the variable Ai ow _k- By fitting cd(Ai ow _]J and C£i(Ai ow _k) to the binding energies 
of 3 H and 4 He, one obtains [56] the values shown in Table 1 for two different choices of the 
low-momentum cutoff. 
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Figure 1: The leading-order contributions to the chiral three-nucleon force: (a) the long-range 
27r-exchange force (b) the medium-range l7r-exchange force and (c) the short-range 

contact interaction vffl- 



Alow-k 


CD 


ce 


2.1 fm- 1 


— 2.062 


-0.625 


2.3 fm" 1 


— 2.785 


-0.822 



Table 1: The values of the low-energy constants cd and ce of the chiral three-nucleon interaction 
fit to the binding energies of A = 3, 4 nuclei for different values of the momentum cutoff Ai ow _ k . 



We now show how to derive from the leading-order chiral three-nucleon interaction, eqs. (1- 
4), an effective density-dependent in- medium NN interaction. We keep our discussion brief and 
refer the interested reader to the original references [44, 57]. We consider the on-shell scattering 
of two nucleons in isospin-symmetric (spin-saturated) nuclear matter of density n = 2k 3 j-/3n 2 in 
the center-of-mass frame, N\{p) + N 2 (— p) — > N% (p + q ) + N 2 (—p — q), where kf is the Fermi 
momentum. The magnitude of the in- and out-going nucleon momenta is \p\ = p = \p + q\ , and 
q — \q\ is the magnitude of the momentum transfer. The resulting interaction in momentum- 
space has the following (general) form 

y{Vi O) — Vc + Ti • r 2 Wc + [Vs + T\ ■ f 2 Ws] ox ■ (?2 + [Vt + ti • r 2 W T ] <?i ■ q o 2 ■ q 
+ [Vso +Tl ■ f 2 W SO ] J r0 2 )-{q*p) 

+ [Vq + Ti-t 2 Wq] ai- (qxp)a 2 - (qXp). (5) 

The subscripts refer to the central (C), spin-spin (S), tensor (T), spin-orbit (SO), and quadratic 
spin orbit (Q) components of the NN interaction, each of which occurs in an isoscalar (V) and 
an isovector (W) version. 

In Figs. 2 and 3 we show the six topologically distinct contributions to the in-medium 
density- dependent nucleon-nucleon interaction to one-loop order. The long-range two-pion ex- 
change three-nucleon force generates the diagrams shown in Fig. 2. Tadpole diagrams obtained 
by contracting the in- and outgoing lines at a tcNN vertex vanish, which leaves only three 
topologically distinct nonvanishing contributions. The short double-line on a nucleon prop- 
agator symbolizes the filled Fermi sea of nucleons, which introduces the "medium insertion" 
— 27i5(k ) 9(kf — \k |) in the in-medium nucleon propagator. Diagram (1) in Fig. 2 represents 
one-pion exchange with a Pauli-blocked pion self-energy correction. The intermediate nucleon 
momentum is independent of the momentum transfer q carried by the pions and therefore 
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Figure 2: In- medium NN interaction generated by the two-pion exchange component of the chi- 
ral three-nucleon interaction. The short double-line symbolizes the filled Fermi sea of nucleons. 
Reflected diagrams of (2) and (3) are not shown. 



the integration over the filled Fermi sea of nucleons produces a correction proportional to the 
density n: 



V, 



med,l 9a M NU 



NN 



n ■ t 2 



a x -qa 2 -q 



(2 Cl ml + c 3 q 2 ) . 



*(ml + qiy 

Since c\ 3 < 0, this term gives a (large) enhancement of the bare l7r-exchange: 



(Itt) _ 9a m n _ _ Q\ ■qo 2 -g 
±611 J* ml + q l 



(6) 



(7) 



In part this enhancement can be interpreted in terms of the reduced in-medium pion decay 
constant, f* 2 s = f 2 + 2c 3 p, which is associated with the spatial components of the axial current. 
Diagram (2) of Fig. 2 represents a vertex correction to one-pion exchange. Fermi sphere inte- 
grals, j|g| <fcF d 3 k, over a single static pion propagator give rise to p-dependent functions Tj(p), 
where the dependence on kp is implicit. For the sake of brevity and continuity in the cur- 
rent discussion, we have left the explicit expressions for the Tj{p) in the Appendix. Summing 
diagram (2) of Fig. 6 together with the three reflected diagrams yields the following contribution 



t rmed,2 
NN 



7T1 • r 2 



(J\-q(J2-q 
m 2 _|_ g2 



O0 , n -> ■ - > > 1 ~ Ac i m l r o(p) + T x (p) 

o2n 6 j^ + q l I L 

-(c 3 + c 4 ) [q 2 (T (p) + 2T 1 (p) + r 3 (p)) + 4T 2 (p)] + 4c 4 



2kj 



- mlT {p) 



(8) 



Kit) 



By analyzing the momentum- and density- dependent factors in eq. (8) relative to V^f N ' , one 
finds that this contribution corresponds to a (large) reduction of the l7r-exchange in the nuclear 
medium. Approximately, this feature can be interpreted in terms of a quenched nucleon axial- 
vector constant g* A . 

Diagram (3) in Fig. 2 represents Pauli-blocked 27r-exchange. The Fermi sphere integrals 
over two pion propagators give rise to p- and g-dependent functions Gj (p, q) which we tabulate 
in the Appendix. Diagram (3) taken together with the corresponding reflected diagram then 
produces the contribution 
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NN 
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-C3 
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+4c 4 fi • t 2 (B x ■ a 2 q 2 - a x ■ q a 2 ■ q)G 2 (p, q) 
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I b Nb 

(4) /%) 




Figure 3: In-medium NN interaction generated by the one-pion exchange component and the 
short-range component of the chiral three-nucleon interaction. 



-(3c 3 + C4T1 ■ f 2 ) i{a x + ff 2 ) ■ (qxp) [2r (p) + 2T 1 (p) - (2m\ + g 2 ) 
x (G?o(p, ?) + 2G?!(p, g))] - Ylc x mli{a x + a 2 ) -(gxp) [g (p, g) + 2G x (p, g) 
+4c 4 f x ■f 2 a 1 -(qx p) a 2 ■ (q x p) [G (p, g) + 4Gi(p, g) + 4G 3 (p, g) 



(9) 



In comparison to the analogous 27r-exchange interaction in vacuum (see section 4.2 in ref. [58]) 
the Pauli blocking in the nuclear medium has generated additional spin-orbit terms, i(a\ + a 2 ) ■ 
(qxp), and quadratic spin-orbit terms, a x ■ (q x p ) a 2 ■ (q x p ), written in the last three lines 
of eq. (9). 

Next we consider the l7r-exchange component of the chiral three-nucleon interaction pro- 
portional to the parameter cd/A x . Contracting two external lines results in two (nonvanishing) 
topologies. By closing a nucleon line at the contact vertex, one obtains a vertex correction to 
l7r-exchange: 



V, 



med,4 



NN 



g A M N c D n _ _ a t -qa 2 -q 
ti ■ r 2 



327r/^A x ' " ml + q* ' (10) 

We note that this contribution is proportional to the nuclear density n due to the g-independent 
Fermi sphere integral. Since Co is negative, V^ d ' 4 reduces the bare l7r-exchange, roughly by 
about 16% at normal nuclear matter density po = 0.16 fm" 3 . Diagram (5) in Fig. 3 represents 
Pauli-blocked (pionic) vertex corrections to the short-range NN interaction. The corresponding 
contribution to the density dependent in-medium NN interaction reads: 



med,5 
NN 



g A M N c D 
64tt 3 /4A x 



Ti ■ T 2 



1a x ■ 02 r 2 (p) + <?i ■ <? 2 2p 2 



— 0^1 • (g x v)°2 



x p 



)) (r (p) + 2r 1 ( P ) + r 3 (p) 



- ax - qa 2 - q I 1 - — I 
+ Ak 3 f - 6mlT (p)\ (11) 



Due to the special spin and isospin structure of this term, it gives nonvanishing contributions 
only to relative S'-states as well as to the 3 Sx — 3 Dx mixing matrix element. 

Finally, we discuss the short-range component of the chiral 3N interaction, represented by a 
three-nucleon contact- vertex proportional to c^/A x . By closing one nucleon line (see diagram 
(6) in Fig. 3) one obtains the following contribution to the in-medium NN-interaction: 



V, 



med,6 _3M N CeU 

~ 87r#A x 



NN 



(12) 



which simply grows linearly in density n and is independent of spin, isospin and nucleon mo- 
menta. 
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Figure 4: Effect (at normal nuclear matter density, no = 0.16 fm -3 ) of the six different terms 
Vnn on the S'-wave components of the free-space nucleon-nucleon potential Viow-fc for a low- 
momentum cutoff of Ai ow _k = 2.1 fm -1 . 



In Fig. 4 we show the contribution from each of these components of the density-dependent 
NN interaction to the 1 Sq and 3 Si diagonal partial wave matrix elements (see ref. [57] for an 
analysis of higher partial waves). The black solid line gives the matrix elements of V\ ow _k 
obtained from the Idaho N 3 LO chiral NN potential for a momentum cutoff of Ai ow _k = 2.1 
fm -1 . The various curves labeled i = 1, 2, . . . , 6 show the matrix elements of Vi ow -fc + V™^ '\ 
We see that the effects from the long-range two-pion exchange part of the chiral three-nucleon 
force are quite large. However, the pion self-energy correction Vjffi' 1 largely cancels the effects 
of the vertex correction V^f' 2 . Since V^^f' 1 and V^ d ' 2 are simply modifications to one-pion 
exchange, a similar cancellation occurs in all higher partial waves as well. The final component, 
V^ d ' 3 , that comes from the long-range three-body force has a much richer spin and isospin 
structure, and at large distances (small momenta) the overall effect is strongly repulsive. 

The two terms V^ d ' 4 and V^ d ' 5 arise from the medium-range one-pion exchange three- 
body force proportional to the low-energy constant cp. In Fig. 4 we see that the effects are 
much weaker than from the long-range chiral three-nucleon force. Since ' 4 is again just a 
reduction of one-pion exchange, it has nonzero contributions in all partial waves. The second 
term arising from the mid-range three-nucleon force, V^ d ' 5 , is also relatively small and gives 
an additional (equal) attraction in the two relative 5-states, as well as the coupled 3 Si — 3 Di 
mixing matrix element not shown in the figure. The final contribution V^ d ' 6 is a pure contact 
interaction and therefore contributes only to the two relative S-states. The interaction is 
spin and isospin independent and therefore contributes to both of these partial waves equally. 
Indeed, the additional repulsion is approximately twice as large as the attraction coming from 
V^ d ' 5 , so taking these two terms together essentially results in a weaker short-range repulsive 
interaction. 



8 



In Fig. 5 we show the density dependence (for n = 0,n /2, and no) of the diagonal S"-wave 
matrix elements of V\ ow -k + Vnn~- Overall, we find that the effects are repulsive and quite 
large in comparison to the free-space low-momentum NN interaction. The same holds in most 
higher partial waves as well, with the exception that the tensor component of the in-medium 
interaction becomes more attractive with the density [57]. The large S'-wave repulsion in the 
three-body forces increases with the density and thereby provides the mechanism for saturation 
of low-momentum interactions [59]. In Section 4 we study the role that these density-dependent 
effects play in describing the anomalously long beta decay lifetime of 14 C and make a comparison 
with analogous interactions obtained from the Brown-Rho scaling hypothesis. 




Figure 5: Partial wave matrix elements for relative S-states from V\ ow -k + Vnn- The low- 
momentum decimation scale is taken to be Ai ow _k = 2.1 fm _1 . 



3 Brown-Rho scaled nucleon-nucleon interactions 

In this section we summarize Brown-Rho scaling and describe a model approach for implement- 
ing these effects in nuclear interactions at densities close to that of saturated nuclear matter. 
The starting point is the scaling relation for (light) hadron masses derived by Brown and Rho 
[12]: _ 

^*=!U 3 (13) 
y g* A M N m a m p m w f n 

which was obtained by implementing scale invariance in chiral effective Lagrangians following 
the approach of Campbell et al. [60]. In eq. (13) the asterisk refers to an in-medium value of 
the given quantity, f* = f* t is related to the time component of the axial vector current, and 
<3> (to be discussed in more detail shortly) is a function of the nuclear density n. Absent from 
the scaling hypothesis in eq. (13) is the pion mass, which in the chiral limit of massless bare 
quarks would be zero due to the Goldstone boson nature of pions. The nonzero pion mass is 
due rather to the explicit breaking of chiral symmetry by the quark mass term in the QCD 
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Lagrangian and is therefore not as sensitive to the background medium as are the masses of 
other hadrons. This is supported empirically from measurements of energy shifts in deeply 
bound pionic atoms [61, 62]. 

From eq. (13), hadronic masses are predicted to scale with the pion decay constant f n , the 
hadronic observable connected to the spontaneous breaking of chiral symmetry. The connection 
between f n and the fundamental order parameter {qq) of QCD is given by the Gell-Mann- 
Oakes-Renner relation [63, 17] 

flml = -m q {uu + dd), (14) 

where m q is the average (current) mass for up u and down d quarks. Assuming that the pion 
mass is protected by chiral invariance, this relation would produce 

where {qq)* = (Q\qq\Q) and denotes the ground state of the dense medium. In the limit of 
low densities the scalar quark condensate is related to the nuclear density by [18, 64] 

{qq) flml 

where ct^n is the pion-nucleon sigma term that describes how the nucleon mass varies with the 
masses of bare quarks. From a dispersion relation analysis of low-energy ttN scattering, one 
can deduce a wN ~ 45 MeV [65]. Thus, at nuclear matter saturation density one would obtain 
from eqs. (13), (15), and (16) that $(n ) ^ 0.83. In Brown-Rho scaling one therefore assumes 
a linear decrease of hadronic mass with the density: 

m* 1 a nN n 

m 2 f 2 m 2 n 

Such a relationship has also be deduced from in-medium QCD sum rule analyses, which can 
provide important constraints on the form the vector meson spectral functions [13, 28, 66]. 

Brown-Rho scaling was one of the first attempts to study - through changes in the symmetry 
breaking pattern of the ground state - the effects of a hot / dense medium on the properties of an 
embedded hadron. Subsequent work had the broader goal to obtain the full spectral function 
of in- medium hadrons, with special emphasis on the vector mesons: p and oj. Within hadronic 
many-body calculations, self-energy corrections S mc d to the meson propagator D(ou,q) defined 
by 

D(u, q ) = 1 (18) 

u 2 - q - m 2 - S vac (g) - £ mcd (w, q ) 

result in modifications to the spectral function A(u, q), which is proportional to the imaginary 
part of the propagator: 

7T (uj 2 — q — m 2 + Re h(u, q)) 2 + 1m q ) 

The real part of the self-energy then contributes to the effective mass of the particle while the 
imaginary part leads to collisional broadening of the width. In a background medium of cold 
nuclear matter, contributions to S me d arise from the coupling of mesons to nucleon-hole and 
resonance-hole excitations [30, 28, 67, 68] as shown in Fig. 6(a). Calculations are simplified 
by invoking the following general result that holds at low-densities (corresponding to small 
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Fermi momenta) and which relates the in-medium self-energy to the total forward scattering 
amplitude T(q, 00) of a vector meson on a nucleon 



S me d(w, q. 



PN 



d 3 p N 
<k t (2n) 3 2E_ 



-T(q,p N ) 



n 



T(q, uo) 



(20) 



where n is the nucleon density and one assumes that the scattering amplitude varies slowly 
with momentum. The result of such studies was that medium modifications to vector meson 
spectral functions in nuclear matter are dominated by subthreshold resonances, such as the 
.Di3(1520) resonance which produces an enhancement in the low-energy region of the p and 
uo meson spectral functions. Most calculations concluded that the spectral functions of vector 
mesons in cold nuclear matter are largely broadened with little shift of the peak mass. However, 
within an in-medium QCD sum rule approach it has been found [66] that the first moment of 
the p-meson spectral function can decrease (similarly to that in eq. (13)) despite the fact that 
the mass peak changes very little. 
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Figure 6: The contribution of resonance-hole excitations to (a) the in-medium uo meson propaga- 
tor, (b) the in-medium NN interaction resulting from uo meson exchange, and (c) the analogous 
contribution to the three-nucleon force. 



Three-nucleon forces and hadronic medium modifications provide two complementary (and 
partly equivalent) frameworks for understanding and implementing density-dependence in the 
nucleon-nucleon interaction. The role of baryonic resonances are important for each: on the 
one hand they give rise to genuine three-nucleon forces, and on the other hand they provide 
the main mechanism for altering mesonic spectral functions in medium. The connection can 
be visualized in the three diagrams of Fig. 6, which show (a) a process contributing to the 
modification of the cu-meson spectral function in medium, (b) the effect that such a process 
would have on the one-boson-exchange NN interaction due to uo meson exchange, and (c) the 
corresponding three-nucleon force obtained by opening the nucleon hole line in diagram (b). 
At present chiral three-nucleon forces are better constrained experimentally and provide a 
complete basis for including effects of the medium, while the properties of hadrons in medium 
are less constrained but have a more direct connection to questions regarding chiral symmetry 
restoration in dense nuclear systems. In the rest of this section we consider a specific model for 
implementing Brown-Rho scaling medium modifications in one-boson-exchange interactions in 
order to make a quantitative comparison to models of the density- dependent NN interaction 
from chiral three-nucleon forces. 

The scaling relation in eq. (13) should be regarded as only an approximate, mean-field 
relation that could be altered through higher-order corrections. Indeed, it is well-established 
[22, 23] that a naive implementation of eq. (13) would yield far too much attraction due to the 
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decreasing of the effective scalar a meson mass. A more accurate treatment of the intermediate- 
range attraction in the NN potential would be to replace the "ficticious" a meson by a pair 
of correlated pions (in the J n = + channel). This has been done in ref. [69], where a meson 
exchange model for the S'-wave nn interaction was developed and found to be dominated 
by the exchange of the p meson as well as repulsive rnr contact interactions constrained by 
chiral symmetry. With this model for the ttti interaction, the nn — > NN amplitude was then 
constructed, which can be related through crossing relations and a dispersion integral to the 
tcN — > irN amplitude from which one can derive the effective 27r-exchange NN potential. 



n/n 


m ai (MeV) 


&/4* 


m a2 (MeV) 


<4/4tt 


0.0 


640 


4.36 


377 


0.640 


0.5 


538 


3.41 


295 


0.450 


0.75 


471 


2.88 


249 


0.309 


1.0 


395 


2.53 


220 


0.264 


1.25 


365 


2.35 


200 


0.234 



Table 2: Masses and coupling constants for the two zero-width bosons parameterizing medium- 
modified correlated 2n exchange in the Bonn-B potential. [23] 

In this formalism finite density effects in the scalar-isoscalar channel arise through the 
density- depedence of the p meson mass {m p ) and form factor cutoff (A p ) as well as the nn contact 
interactions. Moreover, the in-medium pion propagator is dressed with A-hole excitations. In 
practice, this complicated treatment of correlated 2n exchange was modeled with the exchange 
of two scalar-isoscalar bosons with a form factor cutoff of A a = 3 GeV and masses and coupling 
constants shown in Table 2. Together with the scaling of the vector meson masses and form 
factor cutoffs 2 , namely 

=4 = < = £ = i -o.ia.1, (2i) 

m p m w A n 

the modifications to correlated 2n exchange form the basis of the in-medium Bonn-B potential 
[23]. The potential includes also the exchange of the ir, w, and a mesons as described in ref. [71]. 
Using this model of the medium-modified (MM) Bonn-B potential, the properties of symmetric 
nuclear matter were studied [23], and it was found that medium effects in scalar-isoscalar 
exchange diminished with increasing density, thereby allowing for a reasonable description of 
nuclear saturation in a Dirac-Brueckner-Hartree-Fock calculation. 

In Fig. 7 we show the on-shell matrix elements of the low-momentum MM Bonn-B potential 
in relative S-waves for three different densities: n = {0,no/2 and no}. In comparison to the S'- 
wave matrix elements found for the density-dependent chiral interaction, we see that the main 
difference lies in the low-momentum components of the interaction. Whereas the Pauli-blocked 
two-pion exchange diagram in Fig. 3 gives rise to a large repulsion at low momenta (see Fig. 
5), the dropping of the effective G\ and a 2 masses in the MM Bonn-B potential increases the 
medium- and long-range attraction. Indeed, at nuclear matter saturation density, the effective 
o meson masses have dropped to 395 and 220 MeV (from 640 and 377 MeV respectively), 
which extends the range of the attraction significantly. This enhanced attraction is partially 
mitigated by a simultaneous drop in the effective coupling constants g ai and g C2 . However, 

2 To the extent that the phcnomcnological form factors associated nucleon- meson interaction vertices are 
governed by the nucleon radius, they too will decrease with density due to the proposed increase in the nucleon 
size [70]. 
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Figure 7: Diagonal matrix elements for relative S waves obtained from the low-momentum 
medium-modified Bonn-B potential at three densities: n = {0, no/2, no}. 



at higher momenta we see that the two theories lead to similar modifications in the S'-wave 
interactions. In the density-dependent chiral interaction the short-range contact interaction 
proportional to eg provides most of the additional repulsion, whereas in the MM Bonn-B 
potential an analogous effect is achieved through the dropping of the u meson mass. 

One of the original motivations for studying the beta decay of 14 C within the framework 
of Brown-Rho scaling was the sensitivity of the Gamow- Teller matrix element to the nuclear 
tensor force. In fact, employing a residual interaction consisting of only central and spin-orbit 
forces it is not possible to achieve a vanishing matrix element in a pure two p-hole configura- 
tion [45]. Jancovici and Talmi [46] showed that by including a strong tensor force one could 
construct an interaction which reproduces the lifetime of 14 C as well as the magnetic moment 
and electric quadrupole moment of 14 N, although agreement with the known spectroscopic data 
was unsatisfactory. The most important contributions to the tensor component of the nuclear 
potential arise from ir and p meson exchange, which have opposite signs: 



Vj{r) 



VJ(r) 



Np -> c 



47T 



m^-fi • f 2 S 



12 



1 



1 



+ 



3m r ,r 



(m^r) 3 (m^r) 2 "im^r 



(22) 



where S\2 = 3c?i • f a 2 • f — o\ • o% is the tensor operator and /n-k = gAfnn/^f-K- Since the p 
meson mass is expected to decrease substantially at nuclear matter density while the pion mass 
remains relatively constant, an unambiguous prediction of Brown-Rho scaling is the decreasing 
of the total tensor force at finite density. However, as we show in the next section, the additional 
short-range repulsive force resulting from the decreasing u meson mass also plays an important 
role in suppressing the associated Gamow- Teller strength. 
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4 Beta decay lifetime of C 



The extremely long half-life of 14 C (Ti/% ~ 5730 yr [72]) plays an essential role in archaeological 
dating methods, yet naively one would expect the lifetime to be nearly six orders of magnitude 
smaller. The + ground state of 14 C and the 1 + ground state of 14 N to which it decays satisfy 
the selection for allowed Gamow- Teller transitions, and according to the Wigner spin-isospin 
SU(4) symmetry for light nuclei the associated transition matrix element would be unity How- 
ever, in order to reproduce the experimentally observed lifetime, the Gamow- Teller transition 
matrix element must be approximately 2-10~ 3 , which is expected to result from an accidental 
cancellation among otherwise large terms. This transition therefore provides a sensitive test 
for nuclear interaction models and many-body methods. 

Early studies with phenomenological interactions highlighted the importance of the nuclear 
tensor force, without which it would not be possible to obtain a vanishing Gamow- Teller matrix 
element in a 0p~ 2 model space [45, 46]. Employing realistic NN interaction models did not 
improve the theoretical description significantly [73] , but the use of larger model spaces in no- 
core shell model calculations appeared to generally diminish the Gamow- Teller matrix element 
[74]. Despite the long history of the problem, it was only recently suggested that the beta decay 
lifetime is particularly sensitive to the density-dependence of the nuclear force [25, 44]. The first 
of these studies employed the medium modified Bonn-B potential, described in the previous 
section, while the latter made use of the density-dependent chiral NN interaction constructed 
in Section 2. Indeed, since 14 C and 14 N lie just below the double shell-closure at 16 0, the 
density in the vicinity of the valence nucleons is fairly large. In Fig. 8 we plot twice the charge 
distribution of 14 N obtained from electron scattering experiments [75, 76] fit to the harmonic 
oscillator density distribution 

n(r)oc (^ + b^\e- r2 ' d \ (23) 

where for 14 N, b = 1.29 and d = 1.74 fm. Also plotted in the same figure is the square of the 
radial Op wavefunction which peaks at a radius where the nucleon density is close to 90% that 
of saturated nuclear matter no = 0.16 fm~ 3 . The rest of this section is devoted to a systematic 
comparison between the two models of the in-medium nuclear interaction, which despite several 
qualitative differences in their density dependence, both predict that the beta decay transition 
of 14 C is strongly suppressed by medium effects in the nuclear force. 



0.15- 




0.05- 



Figure 8: Twice the charge distribution of 14 N taken from [75, 76] together with the square of 
the radial Op-shell wavefunctions. 
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Figure 9: Diagrams contributing to the effective interaction V e g in our present calculation. 
Wavy lines represent the density-dependent in-medium nuclear interaction caculated in Sections 
2 and 3. 



Nuclei with mass number A > 1 3 are currently beyond the reach of ab-initio methods, such 
as the no-core shell model [77] and quantum Monte Carlo techniques [78]. We approach this 
problem from shell model perturbation theory using a model space consisting of two Op holes 
in a closed 16 O core. The many-body Schrodinger equation 



E \1> 

J -'n Y ni 



is replaced by the model space equation 



(24) 



(25) 



where H e R — Hq + V c q and H is the sum of the kinetic energy and the harmonic oscillator single 
particle potential with oscillator parameter %uj = 14 MeV. We derive the shell model effective 
interaction V e s following the formalism explained in refs. [79, 80]. The effective interaction is 
constructed from a folded diagram expansion in the so-called Q-box, which consists of all irre- 
ducible valence-linked diagrams (usually truncated to a given order). In the present calculation 
we include hole-hole diagrams in the Q-box up to second-order in perturbation theory, shown 
in Fig. 9. In this figure the wavy line represents the full interaction, either VJ ow _fc -j-V™' 1 



NN 



or 



the low- momentum medium- modified Bonn-B potential, V^^ n _ B . Finally, we need the energy 
splitting between the two Op orbitals, e = e(pi/2) — e(p 3 / 2 ), which is most accurately obtained 
from the experimental excitation energy of the first 3/2" state of 15 N located 6.3 MeV above 
the l/2~ ground state. 

We denote the relevant single-hole states by 



|1> 



0^/2 



|2> 



0p-f 2 ) , |3> = |0p-/ 2 



(26) 



and the two-hole states coupled to good angular momentum J, isospin T, and parity tt are 
denoted by \a/3; J n T), where a, (3 = 2 or 3. After constructing the effective interaction with 
the folded diagram method, we diagonalize the following matrix to obtain the eigenvalues and 
eigenstates with J 71 = 1 + , T = (relevant for 14 N): 



(a/3;l+0|K ff |75;l + 0) 




-e 
-2e 



(27) 
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For the J w = + , T = 1 states, we diagonalize a similar 2x2 matrix. The ground states for 
14 C and 14 N are the lowest J n = + , T = 1 and J 71 = 1 + ,T = states respectively, which we 
denote by 



i>i = a 22; 0+1) + b 33; 0+1 



x 



22; l + 0) + y 23; l + 0) + z 33; l + 



(28) 



From these wavefunctions one can derive the following expression for the reduced Gamow- Teller 
matrix element MqY- 



Mgt = E (^/IK^KMII^) = ^= (-2Vbax + 2V2ay + Aby + 2bz) , 



(29) 



where experimentally Mqt — 2 • 10~ 3 . 

For orientation we first determine the low-lying states of 14 C and 14 N, as well as the Gamow- 
Teller transition matrix element connecting the two ground states, using the free-space low- 
momentum N 3 LO potential at a resolution scale of Ai ow _k = 2.1 fm _1 . We will see shortly that 
there is only a small difference when using the free-space Bonn-B potential. For the two + 
states of 14 C we find 



(30) 



V>i(0+,1) = 0.396 |22;0 + ,l) + 0.918 |33;0 + ,1 
ip 2 (0 + , 1) = 0.918 22; 0+ l) - 0.396 33; 0+, 1 



with energy splitting AE(0 + ) = 12.9 MeV. For the two lowest 1 + states we find 



^i(l + ,0) 

<Mi + ,o) 



0.142 
0.361 



22; 0+, 1 ) - 0.681 23; + , 1) + 0.719 33; + , 1 



22; 0+, 1 



0.663 



23; + , 1 



0.656 



33; + , 1 



(31) 



where the first excited state lies AE(1 + ) = 2.76 MeV above the ground state. The second 
excited J w = 1 + ,T = state lies nearly 16 MeV higher in energy than the ground state and 
will therefore be neglected throughout. With these wavefunctions, the ground state to ground 
state Gamow- Teller transition matrix element is found to be 



M GT = -0.90 , (32) 

which is much too large to describe the known lifetime of 14 C. From the ground state wavefunc- 
tions in eqs. (30) and (31) we see that the first three terms in eq. (29) add coherently while the 
last gives a contribution of the opposite sign, though not enough to cancel the first three. In 
order for medium effects to suppress the transition, strength must therefore be shifted to the 
1 33) component of both ground state wavefunctions. 

We now implement density-dependence in the nuclear potential and study the effect on the 
model space interaction matrix elements. In Tables 3 and 4 we show the shell model matrix 
elements for the low-momentum interactions derived from the in-medium chiral N 3 LO and 
Brown-Rho-scaled (BRS) Bonn-B potentials both in free space and at nuclear matter saturation 
density. Table 3 shows the matrix elements for the J w = + ,T = 1 states, while Table 4 shows 
the matrix elements for the J n = 1 + ,T = states. Overall, medium effects in both potentials 
lead to qualitatively similar changes in the shell model matrix elements, though the density- 
dependence of the in-medium chiral interaction is much stronger. This should not be surprising, 
given the stronger S'-wave repulsion in the in-medium chiral interaction, as compared to the 
BRS Bonn-B potential, shown in Figs. 4 and 7. 
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,p = 0+,T = 1 



n/riQ 


(22|V^ al |22> 


<22|V^ai|33> 


(33|V^ d al |33> 





-3.28 


-3.77 


-0.61 


0.25 


-2.64 


-3.66 


-0.05 


1.0 


-0.12 


-2.88 


1.91 




n/n 


(22|V££|22> 


(22 Vg£\M) 


(33|V^r n |33) 





-3.12 


-3.65 


-0.54 


0.25 


-2.72 


-3.38 


-0.33 


1.0 


-1.09 


-2.21 


0.48 



Table 3: Matrix elements (in units of MeV) between 0p~ 2 states coupled to (J W ,T) = (0 + , 1) 
for the in-medium chiral interaction as well as the Brown- Rho-scaled Bonn-B potential at three 
different densities: n — 0, 0.25n , and n . 



r = l+T = 



n/n 


(22|V^ al |22> 


(23|VS5Sa|22) 


(23|\/ c r al |23) 


(23|V^ al |33> 


(22|\/ c r al |33) 


(33IVSSSJ33) 


0.0 


-1.19 


3.89 


-5.20 


1.36 


1.66 


-1.67 


0.25 


-0.61 


3.74 


-4.81 


1.50 


1.33 


-1.00 


1.0 


1.29 


2.81 


-3.03 


1.92 


0.44 


0.72 




n/n 


(22|\/ B m D c n d n |22) 


(23|^r„|22) 


(23|Vffl23> 


(23 Vj^n 33) 


(22|V££|33> 


(33|^|33) 


0.0 


-1.28 


3.50 


-4.91 


1.23 


1.46 


-1.68 


0.25 


-1.09 


3.27 


-4.56 


1.24 


1.38 


-1.57 


1.0 


-0.31 


2.31 


-3.11 


1.26 


1.04 


-1.09 



Table 4: Matrix elements (in units of MeV) between Op -2 states coupled to (J 71 ,T) = (l + ,0) 
for the in-medium chiral interaction as well as the Brown- Rho-scaled Bonn-B potential at three 
different densities: n — 0, 0.25no, and uq. 



At nuclear matter saturation density, medium effects on both potentials are quite large. 
Nevertheless, we can gain insight by studying the interactions at low densities and employing 
perturbation theory. One finds the following corrections to the ground state wavefunctions: 



^ 1} (o + ,i) = vi 0) (o + ,i)- 

- 0.36(33|V mcd |33) ^ UJ (0 + ,1) 



12.9 

Amed| Q Q\l „/,(0) 



0.36(22|\/ mcd |22) + 0.69(22|V mcd |33) 



(33) 



and 

^ 1) (l + ,0) = ^ 0) (l+,0)- 



2.76 



0.05(22|\/ mcd |22) - 0.34(23|V mcd |22) + 0.45(23|\/ med |23) 
- 0.03(23|\/ mcd |33) + 0.17(22|\/ mcd |33) - 0.47(33 1 V" mcd 1 33>] ^ 0) (1 + , 0). (34) 

A straightforward calculation using the matrix elements shown in Tables 3 and 4 for n = 
0.25n gives for the un-normalized perturbed wavefunctions arising from V^? d al 



4 1} (0 + ,1) = C ; (0 + ,l)-0.01^r(0 + ,l) and 



,(o) 
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In-medium chiral N 3 LO 


In-medium Bonn-B 


n/n 


a 


b 


X 


y 


z 


M GT 


a 


b 


X 


y 


z 


M GT 


0.0 


0.40 


0.92 


0.14 


-0.68 


0.72 


-0.90 


0.40 


0.93 


0.11 


-0.62 


0.78 


-0.70 


0.2 


0.38 


0.92 


0.15 


-0.69 


0.71 


-0.91 


0.35 


0.94 


0.08 


-0.55 


0.83 


-0.49 


0.4 


0.36 


0.93 


0.14 


-0.67 


0.73 


-0.82 


0.32 


0.95 


0.05 


-0.50 


0.86 


-0.32 


0.6 


0.33 


0.94 


0.13 


-0.65 


0.75 


-0.74 


0.29 


0.96 


0.03 


-0.45 


0.89 


-0.17 


0.8 


0.30 


0.96 


0.12 


-0.62 


0.78 


-0.64 


0.25 


0.97 


0.02 


-0.40 


0.92 


-0.04 


1.0 


0.26 


0.97 


0.11 


-0.59 


0.80 


-0.53 


0.21 


0.98 


0.01 


-0.36 


0.93 


0.07 



Table 5: The coefficients of the jj-coupled wavefunctions defined in eq. (28) and the associated 
reduced GT matrix element as a function of the nuclear density n for both the Brown-Rho- 
scaled Bonn-B potential and the in-medium chiral nuclear interaction. 

4 1} (l + ,0) = 4 0) (l+,0) + 0.04^ 0) (l + ,0), (35) 

We see that a positive (22; 0+l|\/ mcd |22; 0+1) and a positive (22; 0+l|\/ mcd |33; 0+1) lead to a 
stronger 1 33; + , 1) component in the perturbed 14 C ground state wavefunction. One therefore 
expects quite generally that additional density-dependent repulsion would increase the 1 33; + l) 
component in the 14 C ground state and therefore help suppress the beta decay of 14 C. Both the 
in-medium chiral nuclear interaction and the BRS Bonn-B potential share this feature. 

For the ground state of 14 N the situation is more complicated due to the larger number 
of terms in the perturbed wavefunction of eq. (34) as well as the more complicated depen- 
dence of the shell model matrix elements on the perturbing potential. The two dominant 
terms are (33; l+0|\/ mcd |33; 1+0) and (23; l+0|\/ mcd |23; 1+0), which have the same sign as the 
perturbing potential in ^-waves but enter with opposite sign in the perturbed wavefunction. 
At one-quarter of nuclear matter density, the sum of all six contributions to V^f al yields 
(33; l + 0|V^ d al |33; 1+0) = 0.67 MeV and (23; l + 0|V^ d al |23; 1+0} = 0.39 MeV. The remaining 
four terms together give a minute contribution to the perturbed wavefunction, and therefore we 
expect at low densities for the |33; 1 + , 0) component in the ground state of 14 N to decrease for 
the in-medium chiral interaction and consequently increase the magnitude of the Gamow- Teller 
transition matrix element. Since different components of the in- medium chiral interaction have 
different density-dependences, this behavior can change at higher densities. Indeed by inspect- 
ing the shell model matrix elements at nuclear matter saturation density we see that the opposite 
effect occurs and the GT matrix element is suppressed rather than enhanced (the perturbing 
potential is quite large and therefore our argument based on perturbation theory should be 
treated cautiously). In contrast, the BRS Bonn-B potential has (33; l + 0|Vj^ d 1 33; l + 0) = 0.11 
MeV and (23; l + 0|V^ d |23; 1+0) = 0.35 MeV, which suggests that even at low densities the 
beta decay should be suppressed. 

We now consider the results of the full calculation, including all density-dependent contri- 
butions to second-order in the shell model effective interaction shown in Fig. 9. In Table 5 we 
show the results of our calculations for the expansion coefficients of the ground state wavefunc- 
tions of 14 C and 14 N, as well as the reduced GT matrix element, as a function of the nuclear 
density up to n = n . The full results agree well with the conclusions we obtained from ar- 
guments based on perturbation theory. The Gamow- Teller transition matrix element decreases 
immediately at low densities in the BRS Bonn-B potential, while only after n ~ 0.25n do 
the density-dependent corrections of the chiral nuclear interaction lead to a suppression of the 
matrix element. 
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The Gamow- Teller transition strength, B(GT), is related to the reduced GT transition 
matrix element by 

B(GT) = {g A f-l—\M GT \ 2 ~ |M GT | 2 , (36) 

where Jj = and we have used the approximation that the in-medium axial vector coupling 
constant g* A ~ 1. In general, the effective Gamow- Teller operator has additional terms in a 
nuclear medium [81, 82]: 

G T,eft = 9laL + g* A a + g PA [Y 2 ,a], (37) 

where L is the orbital angular momentum operator and Yi denotes the rank-2 spherical harmonic 
Yi, m - However, from theoretical calculations [83] of the effective GT operator in- medium, as 
well as beta decay calculations [84] performed with phenomenological shell model effective 
interations, it is known that g^A and gpA are almost negligible and that in light nuclei g\ is 
smaller by 15-20% from its free space value of gA = 1-27 (measured in free neutron beta decay). 
In general, one must calculate both the effective interaction V e g and the effective Gamow- Teller 
operator 0GT,eff • In our calculations we assume a 20% reduction of g A in medium and therefore 
set g* A = 1.0. 




Figure 10: The B{GT) values for transitions from low-lying states of 14 C to the ground state 
of 14 N as a function of the nuclear density. The experimental values are from [85]. Note that 
there are three experimental low-lying 2 + states compared to two theoretical 2 + states in the 
Op -2 shell model configuration. 

In Fig. 10 we plot the B(GT) strengths for transitions from the low- lying states of 14 C to 
the ground state of 14 N for both density-dependent NN interactions. These are compared to 
the experimental values [85] obtained from the charge-exchange reaction 14 N(<i, 2 He) 14 C. From 
Fig. 10 one sees that the most dramatic effect is a reduction in the ground state to ground 
state Gamow- Teller transition strength in both models. In the Brown-Rho-scaled Bonn-B 
potential also the first excited (0 + , 1) state of 14 C is decreased toward the experimental value. 
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Figure 11: The splitting between the 1 + ground state of N and the first excited + state as 
a function of the nuclear density for both the medium-modified Bonn-B potential and the in- 
medium chiral nuclear interaction. Also shown is the experimental energy splitting A£ , +-i+ — 
2.3 MeV. 



Such an effect can be obtained from in-medium chiral interactions, but only if the strength of 
the contact three-body force is increased [44]. In the present version of the in-medium chiral 
nuclear interaction, there appears to be very little effect on the other Gamow- Teller transition 
strengths. 

Besides comparing the Gamow- Teller strengths within the two models, it is useful to study 
the energy splitting between the ground state of 14 N and the first excited (0 + , 1) state (the 
isobaric analogue of the 14 C ground state), which is experimentally observed to lie at 2.3 MeV. 
In Fig. 11 we plot this energy splitting as a function of the nuclear density for the in- medium 
chiral nuclear interaction and the BRS Bonn-B potential. It is clear that free-space nucleon- 
nucleon interactions systematically underpredict this level splitting and that additional density- 
dependent repulsion leads to an enhancement of the splitting. The in-medium chiral nuclear 
interaction is much more repulsive than the BRS Bonn-B potential and therefore leads to a 
larger increase in the splitting which overpredicts the value at higher densities. Nevertheless, we 
find that qualitatively the two models for the density-dependent NN interaction are consistent 
with one another and improve the agreement with experiment. 

In order to better understand the relationship between the in-medium chiral nuclear in- 
teraction and the BRS Bonn-B potential, we first study in detail the six different one-loop 
contributions to V^ d arising from the leading-order chiral three-nucleon force. In Tables 6 and 
7 we show the most important shell model matrix elements for each of these components at 
nuclear matter saturation density no- Again, we let perturbation theory arguments direct our 
attention to the relevant physics. As previously mentioned, in order to suppress the Gamow- 
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Teller transition matrix element, it is necessary to increase the 1 33; J n T) components of the 
ground state wavefunctions for both 14 C and 14 N. From eq. (33) it follows that only Vj 
V^ cd ' 3 , and V^ d ' 6 contribute in the desired way. The remaining three contributions to Vj 



med,l 
NN ) 
med 
NN 



are attractive and decrease the strength of the |33;0 + , 1) component of the 14 C ground state 
wavefunction. The structure of the 14 N perturbed ground state is most strongly affected by the 
(33; l + 0|\/ mcd |33; 1+0) and (23; l+0|\/ mcd |23; 1+0) matrix elements. A secondary role is played 
by the (22; l+0|\/ mcd |23; 1+0) and (22; l+0|\/ mcd |33; 1+0) components. From Table 7 we see 
that out of the three contributions that shift strength to the |33;0 + 1) component of the 14 C 
ground state, only the three-nucleon contact interaction proportional to eg also strongly shifts 
strength to the 1 33; l + 0) component of 14 N: 
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(38) 

A full treatment shows that V^ d ' 3 alone has little effect on the ground state wavefunction of 
14 N and in fact actually decreases the strength in the 1 33; l + 0) component of the wavefunction. 
Therefore, the strong short- distance repulsion generated by the three-nucleon contact interac- 
tion is the dominant term contributing to the suppression the Gamow- Teller transition matrix 
element. 



J^ = Q+ ; T = 1 (n = np) 





(22|VgSSa|22> 


(22|V^ d al |33> 


(33|\/ c m h «? d al |33) 


1 


3.88 


0.80 


3.32 


2 


-3.18 


-0.71 


-2.68 


3 


2.15 


0.32 


1.92 


4 


-0.56 


-0.12 


-0.47 


5 


-0.89 


-0.64 


-0.44 


6 


1.75 


1.24 


0.88 



Table 6: Matrix elements (in units of MeV) between Op -2 states coupled to (J 7 ", T) = (0 + , 1) for 
the six density-dependent contributions to the in-medium chiral nuclear interaction at nuclear 
matter saturation density n . 



In the BRS Bonn-B potential one might expect analogously that the scaling of the w-meson 
mass plays a similar role. We therefore have performed a similar calculation using a BRS Bonn- 
B potential with no w-meson mass scaling. In Fig. 12 we have plotted the B(GT) values of 
the ground state to ground state transition using medium-modified interactions with the short- 
distance repulsion removed (the three-body contact term in the case of the in-medium chiral 
interaction and the dropping u meson mass in the case of the BRS Bonn-B potential). We 
find that in both cases there is in fact a strong enhancement of the Gamow- Teller transition. 
Therefore, despite the fact that these two density-dependent interactions are behaving very 
differently at long distances (as inferred from the low-momentum components of their interac- 
tions in Figs. 4 and 7), they both include a strong short- distance repulsion that dramatically 
suppresses the Gamow- Teller transition matrix element. 
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J^ = 1+ ; T = (w = wp) 



n/n 


(22 


(23|V^ al |22> 


<23|V^ r d al |23> 


(23|V^ al |33> 


(22|V^ r d al |33> 


(33|V^ r d al |33> 


1 


4.25 


-0.65 


1.67 


2.14 


-1.67 


3.28 


2 


-3.15 


0.73 


-1.47 


-1.44 


1.47 


-2.65 


3 


1.36 


-0.76 


1.53 


0.00 


-0.86 


1.70 


4 


-0.54 


0.13 


-0.27 


-0.23 


0.25 


-0.47 


5 


-0.47 


0.56 


-1.02 


0.100 


0.14 


-0.35 


6 


1.04 


-1.09 


1.73 


-0.01 


-0.55 


0.88 



Table 7: Matrix elements (in units of MeV) between Op -2 states coupled to (J 7 ", T) = (1 + , 0) for 
the six density-dependent contributions to the in-medium chiral nuclear interaction at nuclear 
matter saturation density n . 



5 Summary and Conclusions 

We have studied the effect of density-dependent low-momentum nucleon-nucleon interactions 
on the beta decay lifetime of 14 C. The in-medium chiral nuclear interaction was obtained from 
the leading-order three-nucleon force by contracting one pair of external lines and summing 
over the filled Fermi sea of nucleons, whereas the Brown-Rho-scaled Bonn-B potential was 
obtained by decreasing the masses and form factor cutoffs of the vector mesons and treating 
microscopically effective a meson exchange as correlated two-pion exchange. Although both 
density-dependent interactions introduce additional strong repulsion at short distances, they 
differ largely at intermediate and long distances. The Brown-Rho scaled Bonn-B potential 
enhances the attraction at these distances (through a decreasing "a" meson mass), but the 
Pauli-blocked two-pion exchange diagram obtained from the long-range chiral three-nucleon 
force is strongly repulsive in this region. Nevertheless, we find that in both models of the in- 
medium nuclear interaction the ground state to ground state Gamow- Teller transition strength 
is strongly suppressed at high densities. In contrast, the GT strengths from the ground state 
of 14 N to the excited states of 14 C exhibit only a small density dependence. From perturbative 
considerations we have found that, in general, short-range repulsive interactions lead to a strong 
suppression of the ground state to ground state transition. Moreover, we have explicitly removed 
the source of short-range repulsion in both interactions and found that this leads rather to an 
enhancement of the Gamow- Teller strength. We therefore suggest that the common short- 
distance behavior of these two interactions is responsible for their qualitatively similar effects 
on the beta decay lifetime of 14 C. 



The authors of this work represent three generations of former students and postdocs who have 
been immensely inspired by Gerry Brown's physics intuition and ideas. We gratefully dedicate 
this paper to Gerry on the occasion of his 85 th birthday. 



22 




Figure 12: The B(GT) values for the 14 C ground state to 14 N ground state transition as a 
function of the density for the medium-modified Bonn-B potential and the in-medium chiral 
nucleon-nucleon potential. The asterisk denotes that the additional density-dependent short- 
range repulsion has been removed from both potentials: the u meson scaling in the case of the 
MM Bonn-B potential and the short-range contact three-nucleon force proportional to ce in 
the in-medium chiral interaction. 
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6 Appendix: Explicit expressions for the Fermi sphere 

/med 
NN 



integrals in V{ ni 



In this appendix we tabulate the different functions that contribute to the density-dependent 
expressions in eqs. (6)-(12). The functions Tj(p) result from Fermi sphere integrals over a single 
static pion propagator and have the explicit forms: 



r o(p) = k f - ra„ 



kf + p k f — p 
arctan — h arctan — 



+ ml+ f'-"\ n m \ + fl +V l , (39) 
4j9 ml + (kf - p) 2 



r i(P) = ^(^+^V)-ro( P )- T ^K + (^ / +p) 2 ] [ml + (kf-pf} In g^j^g , (40) 

r 2 (p) = f + ^/ - - p 2 )T (p) + ±(ml + k 2 f - p 2 )^) , (41) 
fcj m 2 +fc 2 +p2 m 2 + k) + 3p 2 r 

r 3 (p) = ^ ^ r (p) ^ i\(p) . (42) 

The Gj(p,q) functions result from Fermi sphere integrals over the product of two static pion 
propagators. Angular integrations are performed analytically, and the remaining relevant mo- 
mentum integrals are given by 



q Jo , 4^U«2p JA(n\ 



IA(p) + q 2 k 2 ^JA(p) 

where A(p) = [m 2 , + (k + p) 2 ][m 2 + (k — p) 2 }. The remaining functions Gj(p, q) are related to 
those above through the equations 



r , \ r (p) - (ml +p 2 )G (p,q) -G*(p,q) 

= 4p 2 _ g 2 > ( 44 ) 



n , , 3T 2 (p)+p 2 T 3 (p) - (ml+p 2 )G*(p,q) - G„(p,q) 

- 4p 2 _ g 2 ' ( 45 J 

G 2 (p,q) = (ml+p 2 )G 1 (p,q) + G*(p,q) + G u (p,q), (46) 

r( v |£i(p) - 2K + P^ifo g) ~ 2Gi,(p, g) ~ G»(p, q) 

G 3(P,?) = J • (47) 

In this chain of equations the functions indexed with an asterisk play only an auxiliary role 
for the construction of Gi,2,3(p, ?)■ We note that all functions Gj(p,q) are non-singular at 
q = 2p (corresponding to scattering in backward direction). For notational simplicity, the 
/^-dependence of Tj(p) and Gj(p,q) has been suppressed. 
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